Spectral properties of the planar t — J model 
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The single-particle spectral functions A(k,ui) and self-energies E(k,tj) are calculated within the t — J 
model using the finite-temperature Lanczos method for small systems. A remarkable asymmetry 
between the electron and hole part is found. The hole (photoemission) spectra are overdamped, 
with ImE oc lo in a wide energy range, consistent with the marginal Fermi liquid scenario, and in 
good agreement with experiments on cuprates. In contrast, the quasiparticles in the electron part 
of the spectrum show weak damping. 



PACS numbers: 71.27.+a, 79.60.-i, 71.20.-b 

The normal state of superconducting cuprates in many 
aspects contradicts the phenomenology of the normal 
Fermi liquid (FL). Anomalous frequency and tempera- 
ture dependence of several response functions is generally 
attributed to electronic correlations, yet a proper descrip- 
tion is missing so far. The angle resolved photoemission 
(ARPES) experiments Jl]H probe the one-particle spec- 
tral function A(k, uj). At intermediate doping they reveal 
for a wide class of cuprates a well defined large Fermi 
surface (FS) consistent with the Luttinger theorem and 
similar quasiparticle (QP) dispersion |^|. This seems to 
imply the validity of the concept of the usual metal with 
electronic-like FS. Such simple FL picture is in an appar- 
ent contradiction with magnetic and transport proper- 
ties, e.g. electrical conductivity scales with hole concen- 
tration, closer to the picture of holes moving in the anti- 
ferromagnetic (AFM) background. Moreover, in ARPES 
the FL interpretation is spoiled by the overdamped char- 
acter of QP peaks |m- Although a large background 
makes fits of particular lineshapes non- unique J3|,pj , the 
QP inverse lifetime is found to be of the order of the QP 
energy, i.e. t _1 oc uj for uj > T, leading to the concept 
of the marginal Fermi liquid (MFL) (^] with an anoma- 
lous single-particle and transport relaxation, in contrast 
to r _1 cx uj 2 in the normal FL. 

It is unclear whether above features can be reproduced 
within generic models of strongly correlated systems, 
such as the Hubbard and the t — J model, in particu- 
lar in the most challenging regime of intermediate dop- 
ing. Spectral properties of these 2D models have been so 
far studied mainly via numerical techniques [||, e.g. ex- 
act diagonalization (ED) and Quantum Monte Carlo 
(QMC) M. These studies, as well as some analytical ap- 
proaches [|| , established a reasonable consistency of the 
model QP dispersion with the experimental one, as well 
as the possibility of large FS, but have not been able to 
investigate closer the character of QP, being in the core 
of the anomalous low-energy properties. 

The aim of the present work is to employ the finite- 
temperature Lanczos method Q to calculate A(k, uj) 



within the t — J model. This method has been already 
applied to other dynamic [[ll| and static functions, 
yielding features consistent with the MFL concept and 
experiments on cuprates. Although calculations are still 
done in small systems, by using finite (but small) T > 
smooth enough spectra are obtained not only to deter- 
mine the QP dispersion, but for the first time also the 
spectral lineshapes and corresponding self-energies. 
We study the t — J model ^ 



on the planar square lattice and set J/t = 0.3 to address 
the regime of cuprates. The operators Cj S , Cj s project out 
the states with doubly occupied sites. The spectral prop- 
erties of the model Eq. (|l|) are investigated by calculating 
the retarded Green's function (fj, is chemical potential) 
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G(k, W ) = -i / dt e^ + ^({c ks (t),£l(0)}}- 
Jo 



(2) 



The average is grandcanonical, which in actual calcula- 
tions at low T in a system with N sites and fixed hole 
concentration = N^/N is replaced by a canonical one 
in the subspace of states with holes. The two an- 
ticommutator terms correspond at low T to the hole - 
inverse photoemission spectra (IPES), and the electron - 
photoemission spectra (PES), respectively. 

The calculation of G(k, uj) at T = with the ED tech- 
nique is well established |@J^], but a small number of 
sharp peaks in the spectra makes it difficult to extract 
information on lineshapes and self energies. The QMC 
methods resort to the use of maximum entropy analysis 
||, which also leads to quite restricted w-resolution. The 
T > Lanczos method |T(| eliminates these problems 
for dynamic quantities, i.e. yields smoother spectra and 
allows for study of the T-dependence. The requirement 
is however that T > Tf s , where Tf s is the characteristic 
temperature at which in a given small system the finite- 
size effects set in (for discussion of the method we refer to 
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previous works fllO] , |ll[ ). We have calculated the Green's 
function Eq. (0) on systems with N=16 and 18 sites using 
~ 120 Lanczos steps and sampling over ~ 1000 random 
states. The finite-size effects are small at T > T[ S (N, Nh), 
where, e.g., T is ~ 0.1* for N h /N = 3/16. 

From the Green's function we obtain the spectral func- 
tion A(k,(j) = — (l/7r)ImG(k, to) and the one-particle 
density of states (DOS) Af(e) = (2/N) £ k A(k, e - li). 
The latter is used to define the zero of energy and thus 
the chemical potential in Eq. (0) via f^° oo j\f(uo + [i)(e l3uJ + 
l)~ 1 dtu = 1 — Cft,. We find a very good agreement between 
fx calculated this way and from the thermodynamic func- 
tion c h (T,Li) = N h /N g§. 

Of particular interest is the self-energy 



E(k,w) =uj-G{k,ujy 



(3) 



The relation contains no free term, in contrast to the 
usual definition, since the * — J model does not allow for 
a free-fermion propagation even at J = 0. It is also im- 
portant to note that due to projected fermion operators 
in the model the spectral function A(k, uo) is not normal- 
ized to unity (?]], but rather to ({c ks ,c ks }) = (1 + Cft)/2. 
This has several consequences, e.g. ReE(k, uo — » oo) does 
not vanish, but varies linearly with uo. 

In Fig. [I] we first present A(k, uo) for systems with 
Ch ~ 0.12 (combining results for Nh = 2 on systems 
with N — 16,18) and Ch = 3/16. The spectra are 
broadened to Lorentzians of variable width 5 = 6q + 
(<5oo - <5 )tanh 2 (w/A), with S x = 0.2*, 6 a = 0.04*, 
and A = 1.0*. In this way sharper (well resolved) 
low-energy features remain unaffected, while the fluctua- 
tions at higher u>, mainly due to restricted sampling, are 
smoothened out. In any case, S is always smaller than 
the energy scale of main spectral features. 

We observe in Fig. [jj presented at all available k, a 
coexistence of sharper features, associated with coherent 
QP peaks, and of a pronounced incoherent background, 
as already established in earlier studies [0 . The coherent 
peaks in Fig. |l| disperse through uo = as k crosses the 
FS. Within the given resolution in the k-space the FS 
appears to be large already for Ch — 2/18, consistent 
with the Luttinger theorem. The total QP dispersion W 
is broadened as Ch is increased, qualitatively consistent 
with the slave boson picture where W oc c^f + \ J 

In Fig. |^ we show E(k, u>) at c/j = 3/16 and at lowest 
T = 0.1* ~ Tf s . We first notice an asymmetry between 
the PES (w < 0) and IPES {uo > 0) spectra at all k. 
ImE are small for uo > 0, as compared to uo < 0. For k 
outside FS this implies a weak QP damping, consistent 
with sharp IPES peaks seen in A(k,uo), Fig. []} contain- 
ing the major part of the spectral weight. ReE shows 
an analogous asymmetry, in the region uo > resem- 
bling moderately renormalizcd QP. Due to projections in 
Eq. (|l|), the slope in ReE is not zero even at |o;| 3> *, J- 

The behavior on the PES (uo < 0) side is very different. 



For all k, ImE are very large (several * away from u> ~ 0), 
leading to overdamped QP structures. We should here 
distinguish two cases. For k well outside FS, ImE > * 
does not invalidate a well defined QP (at uo > 0), but 
rather induces a weaker reflection (shadow) of the peak 
at uo < 0, as well seen in Fig. |l| for k = (tt, 7t). On 
the other hand, the uo variation for k inside or near the 
FS is more regular, and can be directly related to the 
QP damping. Particularly remarkable feature, found in 
Fig. is a linear frequency dependence of ImE at uo < 
for k = (tt/2, 0), (tt/2, tt/2). Meanwhile k = (0, 0), being 
further away from the FS, seems to follow a different 
(more FL-type) behavior. Such general behavior remains 
similar also for the lower doping c/j = 2/18. 

To address the latter point in more detail, we show in 
Fig. |H the T- variation of ImE for both dopings at selected 
k below the FS. For = 3/16 the linearity of ImE(w) 
is seen in a broad range -2* < uo < at the lowest T 
shown. Moreover, for this higher ('optimum') doping the 
T-dependence is close to a linear one, assuming a small 
residual (finite-size) damping rjo at uo = 0. Data can be 
well described by ImE = ijq + j(\uo\ + £T), with 7 ~ 1.4 
and £ ~ 3.5, baring a similarity to the MFL ansatz [EJ, 
as well as to the conductivity relaxation r^ 1 found in the 
* — J model 1 11 . In contrast, the T-dependence for Ch = 



2/18 seems somewhat different, and ImE oc uo only in the 
interval — * < uo < T. This would indicate the consistency 
with the alternative MFL form , however we should be 
aware that in this 'underdoped' regime finite-size effects 
are larger at fixed T . 

Here we should comment on the manifestation of the 
FS in small correlated systems. At T, uo ~ we are deal- 
ing in the evaluation of Eq.(|^) with the transition be- 
tween ground states of systems with Nh and N' h = Nh ± 1 
holes, respectively. Since these states have definite mo- 
menta ko, they induce strong QP peaks for particular 
k = k^ — ko (defining in this way for a small system the 
FS, apparently satisfying the Luttinger theorem), with 
ImE(k, uo ~ 0) ~ 0. However, the calculated T-variation 
is for a given system meaningful only at T > Tf s . 

From E(k, uo) we can calculate QP parameters: the 
dispersion E^, the weight and the damping Tk, 



E k = ReE(k,£ k ), 

Z k = [l-dIteZ(\L,oo)/doj}zi Ek , 

r u = Z k |ImE(k,£ k )l> 



(4) 
(5) 
(6) 



which are listed in Table 1. We note that parameters are 
of a limited meaning for k inside FS due to large T. In 
particular, £7 k (as well as Z k and r k ) for k = (0, 0) do 
not correspond to a weak QP peak at uo ~ — *, being over- 
whelmed by the incoherent background. Otherwise, the 
enhancement of the dispersion with Ch is seen, accompa- 
nied by a decrease of T for |k| > kp. To establish the re- 
lation with the FL theory one has to evaluate QP param- 
eters at the FS, k = kp. Of particular importance is the 
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renormalization factor Z = Z^ F . Z is still decreasing as 
T is lowered. Nevertheless we find a weak variation (cca. 
20%) within the interval 0.1 < T/t < 0.3, not inconsis- 
tent with the MFL form, leading to Z~ x - \n(uj c /T). 
Regarding the size of Z (at low but finite T > 0) we 
note, that the value of the momentum distribution func- 
tion fits is very close to the maximum for the t — J model, 
n ks - (1 + c h )/2, for all |k| < kp §. Taking the FS vol- 
ume according to Luttinger theorem and assuming that 
nks falls monotonously with |k|, this implies the discon- 
tinuity Z = Sn^g < 2c/j/(l + Ch). We indeed find a con- 
sistent result Z = 0.28 for ch = 3/16, while for c/j = 2/18 
the value is still larger, possibly due to too high T. 

An analogous argument can be used to explain the 
electron-hole asymmetry of A(k, u). Holes added to the 
system at |k| < kp,Lj < move in an extremely corre- 
lated system, strongly coupled to the spin dynamics Q , 
also following the anomalous Iow-lj behavior |fj"l|| . On 
the other hand, states for |k| > kp are not fully popu- 
lated, allowing for a moderately damped motion of added 
electrons for ui > 0. 

Another feature is seen predominantly at smaller dop- 
ing Ch ~ 0.12 for |k| > kp: along with the principal peak 
at lj > a weak bump in the lj < part of the spec- 
trum appears when the FS is crossed along TM. In ReS 
for k = (tt,tt) it emerges even as a strong oscillation, 
leading to a double solution in Eq. (g) In ARPES 

this should be seen as the reappearance of the 'shadow' 
QP band for k above the FS, in accordance with exper- 
iments pj] and some previous studies p5|j7|. The effect 
is less pronounced at larger doping c/j = 3/16, probably 
due to the reduction of the AFM correlation length. 

Finally we show in Fig. |] the variation of the DOS 
AT(s) with doping. For a hole injected in the weakly 
doped system (ch ~ 0.06), a QP coherent peak (of width 



to well understood incoherent hole motion) is dominating 
lower e. At such low doping the electron part of DOS is 
weaker, with the total intensity 2ch as compared to 1 — c/, 
of the hole part. With increasing Ch the hole background 
doesn't reduce in intensity, while the coherent peak near 
the Fermi energy widens and its spectral weight reduces, 
reflecting the broadening of the QP dispersion. At the 
same time, the electron part of DOS is increasing, both 
in the weight and in the width. Note that oscillations for 
e > fi appear in this regime due to the underdamped QP 
and a restricted number of finite-size k. 

Here we mention the relation with the entropy s p2[ , 
assuming the low-T form as follows from the FL theory 
0, i.e. s = ir 2 TJ\f{n)/3Z. With Af{ f i) in Fig. g weakly 
doping dependent at intermediate Ch (also quite close to 
the free-fermion value), and Z ~ 0.28 for Ch = 3/16, we 
get s ~ 0.29 fc^/site at T — O.lt, consistent with static 
calculations [12]. Nevertheless, one should keep in mind 
that such s represents a large increase over the undoped 
(AFM) system, taking into account very few mobile holes 



introduced into the system. 

In the end we comment on the relevance of our results 
to the understanding of the ARPES spectra in cuprates 
0. For lj < we notice the importance of the inco- 
herent background, consistent with the observation that 
in fitting the experiments to either FL or MFL form an 
anomalously large background must be assumed [Q . For 
|k| < kp we find the linewidth typically r ~ t (see Ta- 
ble 0), well compatible (t ~ 0.4 eV in cuprates @|3|) 
with experiments at k away from the FS and at interme- 
diate doping p],^]. Also the MFL form has been claimed 
to describe better the experiments |jj , although this point 
is not yet clarified [|J. We note also that our QP disper- 
sion and the shape of the FS are not entirely of the form 
found experimentally jl| . This could be possibly reme- 
died by including the n.n.n. (f) hopping term [Lj. Still 
we do not expect such corrections to modify conclusions 
concerning the spectral shapes and the QP character. 
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FIG. 1. Spectral functions for Nh = 2 holes on systems 
with N — 16, 18 sites, and for N h = 3 holes on N = 16 sites. 



FIG. 2. Self energy for Ch = 3/16 at various k. 



FIG. 3. ImS(k, lo) for various T for the systems with 
Ch = 2/18 and Ch = 3/16 at selected k below the FS. 

FIG. 4. Single-particle DOS M{e) for variously doped 
systems. For c h ~ 0.06 and c h ~ 0.12 we present joined 
densities for N = 16, 18 systems with Nh = 1 and Nh = 2, 
respectively. The thin vertical lines denote the Fermi energy 
ft. 



TABLE I. QP parameters for two hole concentrations. 

c h = 2/18, T/t = 0.15 c h = 3/16, T/t = 0.15 

k E k /t z k r k /t k E k /t z k r k /t 

(0,0) -3.8 0.80 1.9 (0,0) -4.2 0.73 1.4 

(tt/3, tt/3) -0.7 0.26 0.65 (ir/2,0) -1.1 0.68 1.7 

(2tt/3,0) -0.4 0.35 0.51 (7t/2,7t/2) 0.0 0.28 0.32 

(2tt/3, 2tt/3) 0.5 0.35 0.49 (tt,0) 0.0 0.28 0.32 

(tt,7t/3) 0.1 0.26 0.40 (tt,7t/2) 0.8 0.44 0.31 

(tt.tt) 1.1 0.37 0.54 (tt.tt) 1.7 0.46 0.35 
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